Wave Packet Approach to Neutrino Oscillations with Matter Effects 
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The complete description of neutrino oscillations require the wave packet treatment. For terrestrial 

experiments the contribution from the interaction of neutrino with the Earth matter would modify 

oscillation pattern, i.e., the dependence of the flavor transition probability on baseline L and energy 

E. We thus combine the wave packet approach and the matter effects, in order to have more 

, accurate descriptions of neutrino oscillation. The general expression for the transition probability 

t-H . 

of oscillations with arbitrary numbers of neutrinos are derived. As an example the two-neutrino 
oscillation is studied in detail. 
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I. INTRODUCTION 



iu- 

Over the past few decades the existence of neutrino oscillations has been confirmed in a number of experiments. 
It is widely accepted that neutrinos are massive and mixing. Neutrino oscillations are governed by the mass square 
differences Am^ of neutrino mass eigenstates and the neutrino mixing matrix U a j proposed by Pontecorvo, Maki, 
. Nakawaga and Sakata (PMNS) The standard expression for the probability of the flavor transition v a — > vp in 

qq ■ vacuum is 

o _ _ 

\6 : P ^(L) = \ u ^?Ppj? + 2ReJ2u ai U* j U$ i Up j exp[-i(^ - E. t )L] 

• j j>i 

C\| ' A m 2 L 

~ \U a] ?\U M \ 2 + 2ReJ2u ai U* aj U* 0i Uf 3j exp[-i— ^-], (1) 
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where E and L are, respectively, the average energy of the neutrino beam and the location of neutrino detector with 
respect to the neutrino source. The LjE dependence here is the most important signature of neutrino oscillations 
which have been observed in solar, atmospheric, accelerator and reactor neutrino oscillation experiments. However, 
Eq. (TT]) is a plane- wave approximation and the complete derivation would require the use of the wave packet formalism 
for the evolution of the massive neutrino states. As discussed in refs. 0-3], neutrino oscillations are observable only if 
the process of neutrino production and detection have momentum uncertainties that satisfy the condition of coherent 
production and detection of different neutrino mass eigenstates. Otherwise, different neutrino mass components get 
decoherent and neutrino oscillations are destroyed. Therefore from the quantum mechanical uncertainty principle, 
neutrino states are naturally described by wave packets rather than the plane wave with definite momentum. It is 
revealed that the effects due to such wave packet treatment are observable in oscillations with the existence of sterile 
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neutrinos 6] or oscillation with sufficiently long baseline Q]. In such cases decoherence of different neutrino mass 
eigenstates modify the energy and distance dependence of the flavor transition probability, as one can find in P, Q] 

P aP (L) = £ \U QJ \ 2 \U^\ 2 + 2ReJ2U ai U^U* 0i U pj exp[-2™-i^ - {^f - 2^ 2 (-^) 2 ], (2) 

where L°f c = is the oscillation length, L^° h = 4 ^^ CTx is the coherence length and a 2 = a 2 p + a 2 D is the 

quadratic sum of the production localization and the detection localization. Compared with the standard expression 
Eq. {]]), the two damped terms, — (-rlbr) 2 and — 27r 2 (-^fe) 2 in the exponential, can be understood, respectively, as the 
decoherence due to the separation of different mass eigenstates during propagation, and as the decoherence due to the 
discrimination between different mass eigenstates when neutrinos are emitted and/or absorbed. As any of these two 
terms increases, the oscillation gets suppressed and in the limit, Eq. ((2|) reduces to the averaged transition probability 
Pap{L) = J2j \U a j\ 2 \Upj\ 2 with no more L/E dependence in it. 

In terrestrial oscillation experiments, neutrino beams travel through the mantle and/or the core of the earth. The 
interaction of neutrinos with the particles of matter can change the pattern of oscillations, since the Hamiltonian of 
the neutrinos in matter is different from vacuum. In the flavor basis (y e , i/^, ■ ■ ■ ) T , the effective Hamiltonian is given 
by 

H = U<Ha l g(E 1 ,E 2 , ■■■)U f + diag(l/, 0, 0, ■ ■ ■ ), (3) 

with Ei = \/mf + \p\ 2 is energy of neutrino mass state \vi{p)) in vacuum and V = \plGfn e [n e is the electron 
density in matter) is the charged-current contribution to the matter- induced effective potential of v e 0, [f^] • Note 
that the neutral-current interactions are disregarded here because they are uniform in the flavor basis thus do not 
affect oscillations. Denoting the mixing matrix and Hamiltonian eigenvalues in matter with U and Ej respectively, 
we have 

WHU = diag(£ 1 ,£ 2 ,---). (4) 
Straightforwardly, one could find that the transition probability in matter is simply given by 

P al3 (L) = P a[j {L 1 U -> U, Ej^Ej), (5) 

in the plane-wave approximation. A natural question arises here: what about the complete wave packet approach 
with the existence of matter effects? We will look into it in the following sections. 

II. OSCILLATIONS WITH ARBITRARY NUMBER OF NEUTRINOS 
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It is well known that neutrinos are produced and detected in flavor eigenstates. Considering a production process 

P!^P F + l+ + v a , (6) 
the final neutrino state is generally described with wave packets as 

K) = / d 3 pf P (p,P P )\v a (p)), (7) 
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in which fp(p, Pp) is the momentum distribution function and \v a (p)) is the neutrino state with definite momentum 
p. Usually fp(p, Pp) takes the gaussian form 
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Jp(p,Pp)(xexp[- (P 4 ^ p) ] (8) 

which are sharply peaked around the average momentum Pp with uncertainty up. For process described by Eq. ([6]) 
we have Q 

P P = P Pl - P Pf - P l+ , up = <j 2 Pi + ap F + a? + . (9) 

The first equation is the consequence of momentum conservation and the second one tells that the effective momentum 
uncertainty of the production is dominated by the particle with the largest momentum uncertainty. 

Similarly to the production process presented above, at (L, T) with respect to the source, neutrinos are detected as 
wave packets 

\vp) = J d 3 pf D (p,Pn)\MP)) (10) 

with momentum distribution 

(p-Pd) 

in which Pjj and ap> has analogous definition as in Eq. In order to get the flavor transition amplitude A a p(L, T) = 
(yp\v a {L, T)), we need to perform the space-time evolution operator S = exp(—iHT + ip- L) upon the initial neutrino 
state \v a ). Here the Hamiltonian operator H takes the form of Eq. (|3|) such that matter induced potential V is taken 



■a 



f D oc exp[ —2 ], (11) 



into account. In most of terrestrial oscillation experiments, it was shown that 
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13] 



the relatively little changes 



of electron density n e along the trajectories of neutrinos crossing the Earth mantle or the mantle and the core are 
neglected when the oscillation probabilities are calculated, thus the constant density approximation n ™ an ( COIC ) = 
-man(corc) w jth ^™ an ( corc ) defined as mean electron density number in the mantle (core) is sufficiently accurate. 
Therefore the evolution operator S is invariant as neutrinos traveling through the Earth and the neutrinos arriving 



at the detector are described by 

K(L,T)) = J d 3 pfp(p,P P )exp(-iHT + ipt-L)\is a (p)) 

= f d 3 pf P (p, Pp) Kj exp(-iE s T + ip- L^ip)), (12) 

in which \i>j(p)) denotes the eigenstates of Hamiltonian (|3]). Expanding the detected neutrino state \vp) in Eq. (fT0|) 
with the same eigenstates and making use of the normalization condition (vi(k)\vj (jp)) — SijS^ 3 ' (p— k) we get 

A aP {L,T) = {v p \v a {L,T)) 

= Jd 3 pJ d 3 kfp{p,P P )r D (k,P D )Y J KjUp l eM-^E ] T + ip-L)^ 

= J d 3 pf P (p, P P )f* D (p, Pd) U^Upi exp^T + ip- L) (13) 

j 

For the gaussian distributions Eq. ([8j and Eq. ()11|) it is easy to find that the overall momentum distribution still takes 
the gaussian form with mean momentum P = - and momentum uncertainty = -ir -\ — \-. Notice that 
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both P and a is dominated by the process with the smaller momentum uncertainty. This is due to the fact that a set 
of successive physical processes requires an overlap of the wave packets in momentum space of all the processes, thus 
the one with smallest momentum uncertainty determines the location and shape of the overall wave packet. Therefore 
the amplitude reduces to 

A a0 (L, T) cx f rfVexp[- (p ~f )2 ] £ u* aj U P3 exp(-i^T + ip- L). (14) 

Unlike the vacuum mixing matrix elements U a j which are constants, U a j are functions of the neutrino momenta p. 
However, since the distribution function in Eq. (|14p is sharply peaked around P, it is a good approximation to let 
U a j — U a j | -,_p and take them out of the integration over p. Consequently the integration can be performed with a 
saddle point approximation around P leading to 

A afj (L,T) cx Kfi&i *M-iE<jT + iP-L- ^'/f^ }, (15) 

3 x 

in which E® = Ej\^_p and Vj — -§p\g—p are, respectively, the mean energy and mean group velocity of the corre- 
sponding wave packet in matter. The space-time uncertainty a x = ^ describes localization of both the production 
and detection process because of the relation <t x = g 2 xP + a^D with a xP = an d &xD = are space-time un- 
certainties of the two processes respectively. Notice that opposite to momentum uncertainties, the overall space-time 
uncertainty is dominated by the larger one of the localizations. 

In a majority of neutrino oscillation experiments propagation time T is unmeasured. Even when T is accurately 
measured, the most important observable is the events number accumulated during a period of time. Therefore rather 
than the transition probability P a p(L,T) = \A a p(L,T)\ 2 , we are interested in the time average of it: P a p(L) cx 
J dT\A a f}(L, T)| 2 with the normalization condition ^ Q P a ^(L) = 1. After integration over T we get 

2 
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(16) 



Let us have a look at some features of Eq. (ITBl . compared with the probability in vacuum, Eq. ([2]). The first is, as 
expected, the replacement of mixing matrix elements U a j — > U a j. The square root term can be factored out of the 
summation in relativistic limit thus do not contribute after the normalization. Second, the eigenenergies are changed 
from E® to Ej which is also naturally expected. Moreover, all terms in the exponential are multiplied by combinations 
of the group velocities of the wave packets corresponding to the matter eigenstates as a result, the oscillation 
pattern might be modified depending on the magnitude of matter effect and neutrino energy. In the next section we 
will discuss the two-neutrino oscillation as a specific example of our formalism and then, it will be much clearer about 
how such corrections influence the dependence of transition probabilities on L and E. 



III. TWO-NEUTRINO OSCILLATION 



In general the number of massive neutrinos can be larger than 3, for instance, the existence of sterile neutrinos 
is largely discussed (see [l^] for a brief review). However, all compelling data on neutrino oscillations can still be 
described with three light neutrinos. In a large proportion of experiments, two-flavor neutrino oscillations serve as 
good approximations to the three-neutrino description. Therefore next we discuss the flavor transition between v e 
and within the formalism presented in Sec. II. The Hamiltonian ^ for v e -v^ system is given by 



(17) 




with 



U 



cos V sin ( 
— sin 9 cos ( 



(18) 



It is straightforward to diagonalize such a Hamiltonian and, consequently we get the eigenenergies and mixing matrix 
in matter as 



@l = ~ ( -yj2 {El - E°)Vcos29 - 2£°£° + (E° - £°) 2 + V 2 + + E° + V ) , 



^2 = 2' V 2 ~ v cos 20 - 2E " E 2 + ( E i - E 2? + V 2 + E° + E 2 1 + V) , 



U = 



-+i 



(19) 
(20) 

(21) 



with definitions 



\P\ 



C=[-V+ (E° - El) cos 29 + J (El - E%) 2 + 2F cos 2(9 (£° - £ 2 °) + V 2 csc26>, 



D= V 



E% - Ef) cos 29 + J (El - El) 2 + 2V cos 29 (El - E°) + V 2 esc 29. 



(22) 
(23) 

(24) 



One can verify that our result above is equivalent to the redefinitions of mixing angle in matter 9 m and oscillation 
length in matter Im derived in Ref. 0, [it]] . According to the PREM model 15 1 , the mean electron number densities in 
the mantle and the core in the Earth are respectively n™ an = 2.2 cm _3 NA and n c ° IC = 5.4 cm~ 3 NA- So the magnitude 
of the effective potential inside the Earth is V ~ 10~ 13 eV, which is very small compared with the neutrino energy. 
Therefore we neglect terms of 0(V 3 ). As a result, the appearance probability for the flavor transition v e — > is 
obtained by substituting the eigenenergies El, E 2 , and the mixing matrix elements U a j into Eq. (|16p . which in the 
relativistic approximation leads to 



Pe»(L) 



E^fl^l 2 



2Re{U el U; 2 m i U fi2 exp[-27r J? — - (— K ) 2 - 2^( T ^-) 2 + ( 



r osc v T co 
L 21 iv 21 



k J osc ' 
^21 



- mcoh ' 



(25) 



Here we defined the oscillation length in matter 
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which is naturally expected and the matter-coherence length 



r mcoh — 
J 21 — 



2\f2Ea x 
V sin 26 



(27) 



The matter-coherence term in Eq. (|2"5)) . exp ( riT ^ oh ) 2 , will enhance the oscillation significantly when L approaches 
£2i coh ) while the propagation-decoherence term exp[— (j^k) 2 ] suppresses the oscillation when L approaches L^f 1 - 



10" gV 



which tells 



Considering the magnitude of the mass square difference Am^i ~ 10 eV , we find that " h p- 

L 21 

that the matter-coherence term dominates over the propagation-decoherence term when the neutrino mean energy 

^racoh L C ° h 

E is much larger than 100 MeV and vice versa. On the other hand we have -fhr- ~ 10 eV • a x and -fH^ ~ Ea x , 

^21 Ll 2\ 

which means if we are interested in the matter-coherence effect and/or the propagation-decoherence effect near the 
first oscillation maximum, the localization of either production or detection process must be sufficiently small. And 
for small a x , effects due to the localization-decoherence term exp[— 27r 2 (-^^) 2 ] can be totally neglected. 

L 21 

For oscillations of antineutrinos in matter, the transition probability can formally be obtained by replacing V 
with — V (and <5cp with — Sep in general cases) in the corresponding equations, for instance P g p = P eA1 (y — > — V). 
Notice that the matter-coherence term in Eq. (|25|) enhances both P efl and Pgp, , different from the Mikheyev, Smirnov, 
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3. 



Wolfcnstein (or MSW) effect which leads to resonance enhancement either of P e/J or P g p but not of both 
And due to the dependence of U a j and Ej on V we have neither P eM = Pg^j nor P ee = P ss . This is the result of the fact 
that matter in the Earth is not charge symmetric, thus the effective potential V and consequently Hamiltonian (j3|) is 
neither CP- nor CPT- invariant. 

With the global best fit values of 6*12, Am|i in Ref. [ljj and V(n e — n™ an ), the behavior of P e ^(L) against L and 
E are illustrated in Fig. (JTJ. 
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FIG. 1: Transition probability P(v e -¥ f M ) with a x = 10~ n cm, V = 1.65 x 10" 13 eV, sin 2 612 = 0.312, Amjii = 7.58 x 10" 5 eV 2 
and L — 164 km in the left panel while E — 10 MeV in the right panel. The black solid line denotes the result from Eq. (|25|) . 
the blue dotted line denotes the result of the plane-wave approximation and the red dashed line denotes the result omitting 
the matter effect. 



Let us have a brief analysis of the curves presented above. First, the oscillation of plane-wave approximation is 
apparently larger than the other two. Besides the MSW effect, this is also because the localization a x is set to 
be relatively small, so the decoherence effect due to exp[— (y4rr) 2 ] in Eq. (1231) becomes important, and grows as L 

Lin* 
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increases and/or E decreases. There are also horizontal shift of P e ^{V — > 0) from the other two, which is because of 
the difference between L^ and L^x- Moreover, although the oscillation are both suppressed in P efl and P efJ ,(V — > 0), 
the matter-coherence effect due to exp[— ( j^se ) 2 ] hi Eq. (|2~5l) enhances the probability, such that P e/J is slightly above 
Pe^(V^0) in Fig. ©. 

For oscillation between 2/ e and v T it is just to replace the mixing angle #12 and mass-square difference Am^ in 
each terms of Eq. ([25]) with #13 and Am^. With the best- fit value of #13 from the recent published Daya Bay 
experiment [ijj] and lAm^l from the MINOS experiment [20], the appearance probability P eT of v e o v T oscillation 
with normal mass hierarchy Am^ > is illustrated in Fig. @ and the one with inverted hierarchy Am^ < in 
Fig. ©. 
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FIG. 2: Transition probability P(v e -> v T ) wither = 10" 13 cm, V = 1.65xl0~ 13 eV, sin 2 20 i3 = 0.092, Am^ = 2.32xl0" 3 eV 2 
and L — 1227 km in the left panel while E = 1 GeV in the right panel. The black solid line denotes the result from Eq. (|25|) . 
the blue dotted line denotes the result of the plane-wave approximation and the red dashed line denotes the result omitting 
the matter effect. 
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FIG. 3: Transition probability P(f e 



with o x = 1CT 13 cm, V = 1.65 x 10~ 13 eV, sin 2 20 13 = 0.092, Amfsi = -2.32 x 



10 -3 eV 2 and L = 935 km in the left panel while E = 1 GeV in the right panel. The black solid line denotes the result from 
Eq. (|25p . the blue dotted line denotes the result of the plane- wave approximation and the red dashed line denotes the result 
omitting the matter effect. 



Apparently, the MSW effect leads to the large enhancement in the case of normal mass hierarchy while large 
suppression in the case of inverted mass hierarchy. Besides, the mass- coherence effect and the difference between 
and contribute to both vertical and horizontal shift of P eT from P eT (V —> 0). Similarly to the P e(L1 case, 
as L increases and/or E decreases, contribution of propagation decoherence suppresses the oscillation exponentially. 
Notice that here for v e f-> v T , we need a smaller localization, a x ~ 10~ 13 cm, to have apparent decoherence and 
matter-coherence effect. This is because lAm^l and E are both larger than the previous case, and in order to avoid 
L§° h ^ Lfi and/or L^f oh 3> L^\ c , smaller localizations of both production and detection are necessary. 

We conclude this section with some comments on the main results Eq. (1251) . Fig. ([1]), Fig. @ and Fig. ([3]). The 
two neutrino flavor transition probability with matter effects in the wave packet formalism is accurately described by 
Eq. (|25p . as long as the matter induced effective potential, V, is small, since we have ignored contributions of 0(V 3 ). 
Otherwise, the matter-coherence term exp[( rn ^ oh ) 2 ] would destroy the unitarity of probability unless we include higher 
order terms of V. In all the three figures, we present the behaviors of PepM > ^e^(r) (V ~> 0) and Pe^(r) with plane- 
wave approximation, from which we find that if oscillation parameters (L, E and especially a x ) are set properly, the 
predictions of wave packet approach with matter effects are apparently different from the predictions either without 
matter effects or without wave packet effects. Therefore, if observations on the wave packet formalism are performed, 
contributions of matter effects can not be trivially ignored. We note here that the two-neutrino oscillations serve 
only as approximations of the physical truth. For accurate data analysis, it is required to work with three-neutrino 
oscillations which, mathematically, is a little bit more complicated. 

IV. DISCUSSIONS AND CONCLUSIONS 

In the wave packet formalism propagation decoherence is fundamentally due to the difference between the group 
velocities of the different mass components. For the two-neutrino oscillations with matter effects discussed in the 
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previous section, we have 



Vi 



V 2 sin 2 26> 



IE 2 



i' 2 



V 2 sin^ 26» 



2E 2 



(28) 



for velocities with Amfi > 0. It is indicated that matter effects decrease the group velocity of the lighter neutrino 
wave packet while increase the group velocity of the heavier neutrino wave packet. Notice that when neutrino energy 
is sufficiently large, more terms need to be included and the velocities will be always less than 1. 
Recent reactor oscillation experiments 
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21( have revealed that 6*13 is much larger than 0, which makes the 



measurement of CP violation in neutrino oscillations possible and easier. Besides the contribution of the CP phase 
#cp in the mixing matrix, there are also CP violation due to the matter effects, even in the two-neutrino oscillation 
system as discussed in the previous section. As an example, defining CP asymmetry A~F = P eii — Pep,, Fig. (j4]) shows 
the behavior of compared with the plane-wave approximation. 
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FIG. 4: CP asymmetry with a x = KT 11 cm, V = 1.65 x 10~ 13 eV, sin 2 0i 2 = 0.312, Am^ = 7.58 x 10~ 5 eV 2 and 
L — 164 km in the left panel while E — 10 MeV in the right panel. The black solid line denotes the result from Eq. (I25[) and 
the blue dotted line denotes the result of the plane-wave approximation. 



As indicated in Fig. (j4j, when L increases or E decreases the difference between the wave packet result and the 
plane-wave approximation becomes more significant. Therefore, in order to measure CP violation precisely, it is 
necessary to work within the wave packet formalism in certain cases. 

In conclusion, for terrestrial experiments, the wave packet description of neutrino oscillations, Eq. @, is modi- 
fied by the interactions between neutrinos and the Earth matter they travel through, as indicated in Eq. (1161) and 
Eq. (f25|) . which are the main results of this paper. From Fig. (JT]), Fig. ([2]) and Fig. ([3]), we find that the complete wave 
packet approach with matter effects to neutrino oscillations gives distinctive predictions under certain parameter set- 
tings. Therefore, for future terrestrial neutrino oscillation experiments, especially with small localizations of neutrino 
emission and absorption, wave packet approach with matter effects can describe the data with more accuracy. 
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